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B0850+054: a new gravitational lens system from CLASS 
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1 INTRODUCTION 



The Cosmic Lens All-Sky Survey (CLASS) (Myers et al. 2002 



ABSTRACT 

We report the discovery of a new gravitational lens system from the CLASS survey. Radio 
observations with the VLA, the WSRT and MERLIN show that the radio source B0850+054 
is comprised of two compact components with identical spectra, a separation of 0.7 arcsec and 
a flux density ratio of 6:1. VLBA observations at 5 GHz reveal structures that are consistent 
with the gravitational lens hypothesis. The brighter of the two images is resolved into a linear 
string of at least six sub-components whilst the weaker image is radially stretched towards the 
lens galaxy. UKIRT if -band imaging detects an 18.7 mag extended object, but the resolution 
of the observations is not sufficient to resolve the lensed images and the lens galaxy. Mass 
modelling has not been possible with the present data and the acquisition of high-resolution 
optical data is a priority for this system. 
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servations of B0850+054 in the radio, optical and infra-red which 
leave little doubt that B0850+054 is a lens system. 



Browne et al 



2002) is a survey of flat-spectrum sources with a flux 



density in excess of 30 mjy at 5 GHz. 22 lenses have been found 
in CLASS from a total of 16,503 sources observed. The main aim 
of this survey is to identify gravitational lens systems. These can 
be used_^ocojistrau^he_cosmoloe^ 



be userjJocmistnvm^he_cosir^^ 

et al. |199S|; Fassnacht et al||l999t iKoopmans et alj|2000| ) , f2 and 



et al. UWft passnacnt et aljjivy^ iRoopmans et alJlZUUU ) i Z« a 
A (e.g. |Kochanck||l996| ; tlelbig et alj |l999fc IChae et afl|2002 l 



well as galaxy mass profiles (e.g. 
200ll). 



Rusin & Ma 2001 



Cohn et al 



CLASS B0850+054 has a 5-GHz GB6 ( |Gregory et all|l996| ) 
flux density of 68 mjy and a two-point spectral indexfl between 1 .4 
and 5 GHz of a = —0.2, thus satisfying the criteria for inclusion in 
the CLASS sample. The original 8.5-GHz VLA snapshot (30-s on 
source) revealed two compact components separated by 0.7 arcsec, 
with a flux density ratio of ~6:1. Here we present follow-up ob- 
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We define the spectral index, a, such that S 



2 RADIO OBSERVATIONS 

B0850+054 has been observed with the VLA, the Multi-Element 
Radio-Linked Interferometer Network (MERLIN), the Very Long 
Baseline Array (VLBA) and the Westerbork Synthesis Radio Tele- 
scope (WSRT). A summary of the observations presented in this 
paper is given in Table [l| 



2.1 VLA 

The VLA data consist of observations at 5, 8.5 and 15 GHz. At 
each frequency two channels of 50-MHz bandwidth were used and 
the target observed for about 10 min. Initial calibration was done 
with the NRAO AIPS pack age and the data were then mapped and 
self-calibrated in DIFMAP (Shepherd 1997). The resultant maps are 
shown in Fig. |l|. These show that the source consists of two com- 
pact components aligned approximately north-south. As is custom- 



© 0000 RAS 



2 A. D. Biggs et al. 



Table 1. Radio observations of B0850+054 presented in this paper. 



Date 


Array 


v (GHz) 


1998 Mar 15 


VLA 


8.5 


2000 Oct 07 


VLBA 


5.0 


2001 Jan 27 


MERLIN 


5.0 


2002 Feb 04 


WSRT 


0.360 


2002 Apr 14 


VLA 


5.0 


2002 Apr 14 


VLA 


8.5 


2002 Apr 14 


VLA 


15.0 



ary, the brighter of the two components is denoted image A and the 
other image B. The results of model-fitting to the (u, v) data in 
DIFMAP, using an unresolved point (delta) component to represent 
each image, are presented in Table [|. We also tabulate the flux den- 
sity ratios, demonstrating that the spectral index of each component 
is the same. 



2.2 MERLIN 

The MERLIN 5-GHz observations were taken over a period of 9 h, 
with a bandwidth of 16 MHz, and were phase-referenced. The data 
were calibrated in AIPS and mapped and self-calibrated in DIFMAP; 
a map made from naturally-weighted data is shown in Fig. As 
with the VLA data we modelfit using DIFMAP and show the results 
in Table ^ The higher resolution of MERLIN allows a more ac- 
curate estimation of the image separation (679 ± 1 mas) and also 
allows us to detect signs of resolution in image A. For this reason 
we fit an elliptical Gaussian to image A whilst still modelling im- 
age B with a delta component. No evidence for polarisation of the 
images has been found down to the noise limit. This implies that 
the polarisation of the lensed radio source is < 1 per cent (5 a) at 
5 GHz. 



2.3 VLBA 

Phase-referenced observations were taken with the VLBA at a fre- 
quency of 5 GHz for a total of 12 h. A total bandwidth of 32 MHz 
was observed in four separate, but contiguous, bands of 8 MHz 
each. The data rate was 128 Mb s~ . Calibration and mapping 
took place within AIPS following standard procedures. Natural! 
weighted maps of images A and B of 0850+054 are shown in Fig 
Image A is found to have a total length of about 20 mas and con- 
sists of a bright "core" region (containing two sub-components of 
similar brightness) near to its eastern extremity and several bright 
knots to the west of this. To the east of these sub-components there 
is a short extension which is most likely a counterjet. As expected 
for a de-magnified version of A, image B is relatively compact. It 
is, however, resolved north-south suggesting that with higher reso- 
lution the sub-components in the jet could be mapped. 



2.4 WSRT 

In order to obtain a low-frequency flux density measurement of 
B0850+054, observations were made with the WSRT with its 
broadband "92-cm" system for two periods of 2 h, separated by 
about 4 h. Three of the eight 5-MHz bands were affected by RFI 
and discarded. After mapping and self-calibrating we detected a 
source of 50 ± 5 mjy close to the position of the target. Due to the 
very elongated beam at this declination the possibility of confusion 
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Figure 3. UKIRT K-bwd image of B0850+054. Contours are plotted at the 
background sky level plus linearly-increasing multiples (—1, 1, 2, 3, etc) of 
three times the off-source rms noise. The crosses mark the position of the 
A and B components. 



with un related sources is si gnificant. Af ter searching the 1.4-GHz 
NVS S ( |Condon et al1]l998| ) and FIRST ( [Becker, White & Helfand 



1995) surveys, both of which detect B0850+054, we identified a 
~14 mjy confusing source 22 arcsec west and 144 arcsec north of 
B0850+054, well within the WSRT beam. As this source is also de- 
tected in our VLA data at 5 and 8.5 GHz, together with the 1 .4-GHz 
flux densities we are able to calculate its spectral index, a — —0.8. 

Assuming that the spectrum of this source does not turn over, 
its flux density at 360 MHz would be ~41 mjy, thus suggesting that 
it is contributing significantly to the measured WSRT flux density. 
We have aligned the Right Ascension reference frame of the WSRT 
360-MHz image with that of the NVSS using six sources that are 
common to both. The scatter in this relative astrometry is less than 
1 arcsec and from it we have determined that the WSRT detection 
lies approximately midway between the NVSS positions of the two 
sources. Therefore we conclude that each contribute about equally 
to the combined flux density and adopt a value of 25 ± 5 mjy for 
B0850+054 at 360 MHz. 



3 UKIRT OBSERVATIONS 

Service observations with the United Kingdom Infra-Red Tele- 
scope (UKIRT) in K band (2.2 /im) were taken on the nights of 
2001 March 17 and 2001 March 18 using the UKIRT Fast-Track 
Imager (UFTI). Both nights' observations were made using a mo- 
saic mode, each mosaic consisting of nine 60 s pointings with an 
offset between each of 10 arcsec. A total of nine mosaics were ob- 
tained over the two nights. Short observations (20 s) of the stan- 
dard star FS15 were taken on each night so that the data could be 
flux calibrated. Data reduction was performed using the UKIRT 
pipeline analysis tool, ORAC-DR. The J^-band image resulting 
from the sum of all mosaics is shown in Fig. j| 

We detect an 18.7 mag source at the location of the radio im- 
ages, but unfortunately the seeing is not adequate to resolve the 
radio images. From measurements of the sizes of stars in the field 
we measure a seeing-disc size of about 0.6 arcsec, only marginally 
less than the separation between the radio components. The source 
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Table 2. Results of model-fitting to the MERLIN and VLA data. The VLA data were model-fitted assuming delta components for both A and B whilst the 
MERLIN data were fitted with an elliptical Gaussian to A and a delta component to B. Included are the flux densities of images A and B, their flux ratio, 
their separation (r) and their orientation (9) measured in degrees east from north. We adopt errors of 3 per cent on the VLA 5 and 8.5-GHz flux densities and 
5 per cent on the MERLIN 5-GHz and VLA 15-GHz. 



v (GHz) 


Array 


S„,A (mJy) 


S„,b (mJy) 


S'i/.a/'Sv.b 


r (mas) 


en 


5 


MERLIN 


49.0 


8.1 


6.0 ±0.4 


679 ± 1 


166.5 ±0.1 


5 


VLA 


55.1 


9.0 


6.1 ±0.3 


683 ±8 


166.2 ±0.7 


8.5 


VLA 


40.4 


6.6 


6.1 ±0.3 


679 ± 5 


166.6 ± 0.4 


15 


VLA 


26.2 


4.4 


6.0 ±0.4 


680 ±8 


166.6 ±0.7 




Figure 1. VLA and MERLIN maps of B0850+054. Top left: VLA 5-GHz map with a restoring beam of 524 X 356 mas at a position angle of -49?9. The 
data are uniformly weighted and the rms noise in the map is 157 /^Jy beam -1 . Top right: VLA 8. 5-GHz map with a restoring beam of 317 X 200 mas at a 
position angle of —51? 6. The data are uniformly weighted and the rms noise in the map is 121 /iJy beam -1 . Bottom left: VLA 15-GHz map with a restoring 
beam of 153 X 123 mas at a position angle of — 46?2. The data are uniformly weighted and the rms noise in the map is 288 fi]y beam -1 . Bottom right: 
MERLIN 5-GHz map with a restoring beam of 114 X 38 mas at a position angle of 22? 2. The data are naturally weighted and the rms noise in the map is 
100 /iJy beam -1 . All maps are plotted on the same angular scale and each has the restoring beam plotted in the bottom-left corner. Contours are plotted at 
multiples (—1, 1, 2, 4, 8, 16, etc) of 3<r where a is the off-source rms noise in the map. Image A is the northern of the two components. 
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Figure 2. VLBA 5-GHz maps of image A (left) and image B (right) of B0850+054. The maps are uniformly weighted and both are plotted on the same angular 
scale. The restoring beam is shown in the bottom-left of each map and has a FWHM of 3.1 X 1.6 mas at a position angle of — 3?8. Contours are plotted at 
multiples (—1, 1, 2, 4, 8, 16, etc) of 3<r where a is the off-source rms noise in the map (50 /ijy beam^ 1 ). The grey scale represents surface brightness in units 
of mjy beam" 1 . Positions are offset from 08 h 52 m 53f 5772, +5°15'15'.'329 (J2000). 



is however extended, with a major axis (FWHM ~ 0.8-0.9 arc- 
sec) that is aligned approximately north-south {9 ~ 10°). As this 
is similar to the orientation of the lensed images it is likely that 
some lensed emission has been detected. We have attempted to ac- 
curately align the radio and UKIRT images by using the position of 
one of the stars in the UKIRT field that has been measured with the 
Catisberg Automatic Meridian Circle (D. W. Evans, private com- 
munication). The positions of the radio components are also plot- 
ted in Fig. pi The accuracy of our astrometry is set by the positional 
accuracy of the Carlsberg star which, with a positional uncertainty 
of about 150 mas, means that we cannot be certain if most of the 
emission we see is from the lens or the lensed images. 



4 DISCUSSION 

We have presented radio and infra-red observations of the CLASS 
source B0850+054. In this section we describe the evidence that 
leads us to believe that this source is in fact a doubly-imaged lensed 
system. 

The most convincing evidence is provided by the radio data. 
Our VLA, MERLIN and WSRT measured flux densities at 15, 8.5, 

5 and 0.325 GHz, plus the 1.4-GHz NVSS measurement are plotted 
as a function of frequency in Fig. ^. Below 5 GHz the existing maps 
do not have the resolution to enable us to separate the A and B im- 
ages and so only a total A+B flux density is plotted. At frequencies 
^5 GHz the excellent agreement between the radio spectra of A 
and B is clear and the flux ratios in Table [2] confirm that the two 
components have identical spectra to within the errors. 

Further support for the lensing hypothesis comes from a con- 
sideration of surface brightness arguments. As a lens system pre- 
serves surface brightness, the weaker image must be the more com- 
pact as demonstrated particularly well by the VLBA data. More- 
over, the extension of B along the line of the A-B separation is 
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Figure 4. Radio spectrum of B0850+054. The solid line shows total flux 
densities (A+B) whilst below this the individual A and B flux densities are 
plotted where possible. 



characteristic of gravitational lensing and the lack of detectable po- 
larisation in each image is also consistent with the lensing hypoth- 
esis. 

A spectrum taken with the W. M. Keck I telescope has mea- 
sured a redshift for the lensing galaxy of zi = 0.59 (McKean et al., 
in preparation), based on the detection of Ca H & K and G-band 
absorption as well as a 4000 A break. Also detected is a single 
emission line which is almost certainly associated with the lensed 
radio source. The most likely identification for this line is Mg II 
(z s = 1.14) though it could possibly be Lya (z s — 3.93). As- 
suming that Mg II is the correct identification, the image separation 
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together with the measured lens redshift and a source redshift of 
1.14 predicts a K magnitude for the lens galaxy of 17.9,^]close to 
the measured UKIRT magnitude of 18.7. Unfortunately, the present 
optical data do not provide enough constraints (in particular, the 
galaxy position relative to the lensed images is unknown) to make 
any attempt at mass modelling fruitful. 

From Fig. ^ it can be seen that the source spectrum has a turn- 
over at approximately 1 GHz, thus making it a member of the giga- 
hertz peaked spectrum (GPS) class. Low flux density variability is a 
generall y recognised characteristic of GPS sources, as is low polar- 
isation (O'Dea 1998) and so B0850+054 seems consistent with the 



GPS identification. The source is < 1 per cent polarised, and would 
also appear to be only moderately variable, if at all. Flux density 
measurements from the VLA and MERLI N data, the GB6 su rvey 
and the Parkes-MIT-NRAO (PMN) survey ( [Griffith et alT]l995| ), all 
at 5 GHz, are identical to within their errors. At 8.5 GHz, recent 
monitoring with the VLA (six epochs over three months) found no 
evidence for significant variability. Furthermore, an intrinsic size 
for the GPS source of ~10 mas (as measured from the VLB A map 
of image A and corrected for a lens magnification of ~2, a plausible 
value) is consistent with the frequency and flux density (also cor- 
rected fo r the lens magnificat ion) of the turnover using the relation 
found by 



Snellen et al 



(200C) for GPS sources. 



5 SUMMARY AND FUTURE WORK 

We have discovered a new lens system during the course of the 
CLASS survey, B0850+054. It consists of two images of a GPS 
radio source with a separation of 0.7 arcsec and a flux density ratio 
of ~6:1. UKIRT K-band imaging does not resolve the lens galaxy 
and the lensed images, but a lens galaxy redshift of 0.58 has been 
measured with the W. M. Keck telescope. 

A priority for future work is the acquisition of a high- 
resolution optical/infra-red image either with the HST or with 
ground-based adaptive optics. Knowing the relative positions of the 
lens and the lensed images will enable detailed mass modelling to 
be undertaken. We also plan higher resolution VLBI observations 
to look for the substructure in image B that must be there. At the 
same time we hope to detect motion of the sub-components. There 
is the exciting possibility, if B0850+054 is really a two-sided ob- 
ject, of being able to measure the advance speeds of the approach- 
ing and receding jets, made more easily visible by the lens mag- 
nification. When combined with measurement of the relative flux 
densities of the receding and approaching components, this opens 
up th e prospect of an independen t constraint on the Hubble constant 
(e.g. Taylor & Vermeulen 1997). 
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Invoking the Faber-Jackson relation (Faber & Jacksor I97f ) and assum- 



ing that an L* galaxy at z = 0.6 has a K magnitude of 16.4. 
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